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The formation and evolution of the elliptic flow pattern in Pb + Pb collisions at y/s = 5.5A TeV 
and in Au + Au collisions at y/s = 200A GeV are analyzed for different hadron species within the 
framework of the HYDJET++ Monte Carlo model. The model contains both hydrodynamic state 
and jets, thus allowing for a study of the interplay between the soft and hard processes. It is found 
that jets terminate the rise of the elliptic flow with increasing transverse momentum. Since jets are 
more influential at the Large Hadron Collider (LHC) than at the Relativistic Heavy Ion Collider 
(RHIC), the elliptic flow at LHC should be weaker than that at RHIC. The influence of resonance 
decays on particle elliptic flow is also investigated. These final state interactions enhance the 1ow-pt 
part of the i>2 of pions and light baryons and work toward the fulfillment of idealized constituent 
quark scaling. 
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I. INTRODUCTION 

In noncentral collisions between two nuclei the beam 
direction and the impact parameter vector define a re- 
action plane for each event. The observed particle yield 
versus azimuthal angle with respect to the event-by-event 
reaction plane might carry information on the early col- 
lision dynamics [1, 2]. An initial nuclear overlap region 
has an "almond" shape at nonzero impact parameter. If 
the produced matter further interacts and quickly ther- 
malizes, the pressure that is built up within the almond- 
shaped region, develops anisotropic pressure gradients. 
This pressure pushes against the outside vacuum and the 
matter expands collectively. The result is an anisotropic 
azimuthal angle distribution of the detected particles. 
One can expand this azimuthal angle distribution in 
Fourier series [3, 4]. The second coefficient of the ex- 
pansion V2 is called elliptic flow: 

^<cos20) = /4^4\ ■ (1) 

\Px+Py I 

Here, <fi is the azimuthal angle of a particle relative to 
the reaction plane, and p x and p y are the in-plane and 
out-of-plane components of the transverse momentum of 
a particle, respectively. 
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It was found [2, 5] that anisotropic flow is a 
sclf-qucnching phenomenon since it reduces spatial 
anisotropy as it evolves. Therefore, the observed elliptic 
flow must originate at early stages of the collision when 
the anisotropy is still present in the system, and no flow 
is generated when the spherical symmetry is restored. 
Thus, the elliptic flow can provide information about hot 
and dense matter created in relativistic heavy-ion colli- 
sions. 

Although the elliptic flow has been extensively stud- 
ied both theoretically and experimentally (see, e.g., 
Refs. [6, 7] and references therein), its behavior in Pb 
+ Pb collisions at Large Hadron Collider (LHC) energy 
y/s = 5.5A TeV, compared with that at the Relativistic 
Heavy Ion Collider (RHIC), y/s = 200 A GeV, remains 
completely unclear (for a recent review see Ref. [8]). 
Some models predict further increase of the v-i at LHC, 
while others favor a similarity, ss v^ m , or even a 

decrease, v\ nc < v^ mc . The difference among the pre- 
dictions of different models comes from the treatment of 
various processes at partonic and hadronic levels, such as 
equilibration, features of quark-hadron phase transition, 
hadronic cascade, equation of state (EOS), and cross sec- 
tion of partonic scattering, as well as from the initial 
conditions. To the best of our knowledge, the interplay 
between the ideal hydrodynamic behavior and jets, that 
is, the influence of the jets on hydrodynamic flow with 
rising collision energy, has not been elaborated yet. For 
this purpose we employ the HYDJET++ model [9] to 
estimate the azimuthal anisotropy of particles in Pb + 
Pb collisions at y/s = 5. 5 A TeV and compare the ob- 
tained results with calculations for Au + Au collisions at 
the top RHIC energy y/s = 200A GeV and also with the 
experimental data. 

The HYDJET++ model [9] is a superposition of soft 
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and hard parts. These parts are independent and their 
contribution to the total multiplicity production depends 
on collision energy and centrality and is tuned by model 
parameters. The hard part of the model is identical 
to that of the HYDJET model [10] and can account 
for jet quenching effect and the shadowing effect [11]. 
The soft part of the HYDJET++ event represents the 
"thermal" hadronic state where multiplicities are deter- 
mined under assumption of thermal equilibrium [12, 13]. 
Hadrons are produced on the hypersurface represented 
by a parametrization of relativistic hydrodynamics with 
given freeze-out conditions. The model is capable of si- 
multaneously reproducing the main features of heavy-ion 
collisions at RHIC, namely, (i) hadron spectra and ratios, 
(ii) radial and elliptic flow, (iii) femtoscopic momentum 
correlations, as well as (iv) high-py hadron spectra. 

As the fireball expands, its temperature and energy 
density drop. Finally, at the freeze-out stage the system 
breaks up into hadrons and their resonances. The effect 
of resonance decays, namely, final state interactions, on 
the resulting elliptic flow of particles is quite important 
at both LHC and RHIC energies. Here, we are going to 
benefit from the rich table of baryon and meson resonance 
states implemented in HYDJET++. 

The paper is organized as follows. Features of the 
model are presented in Sec. II. In Sec. Ill the influence 
of jets on the elliptic flow of the most abundant hadrons, 
such as pions, protons, kaons and lambdas, is studied. 
We show that the superposition of hydrodynamic flow 
and jets leads to reduction of the flow at high trans- 
verse momenta. Compared with RHIC, the jet fraction 
at LHC starts to dominate at smaller pt, thus effectively 
decreasing the elliptic flow of all particles. In Section IV 
we describe the study of the influence of resonance decays 
on the V2 of stable (with respect to strong interaction) 
particles. At both RHIC and LHC, the effect is found 
to be the strongest for protons and lambdas, relatively 
moderate for pions, and almost absent for kaons. While 
jets modify the (semi)hard part of the V2(pt) spectra, the 
contributions from resonance decays alter their soft part. 
In Sec. V we show that these contributions can account 
for better realization of the constituent quark scaling of 
the hadron elliptic flow. Finally, conclusions are drawn 
in Sec. VI. 



II. THE HYDJET++ GENERATOR 

HYDJET++ is a Monte Carlo event generator for the 
simulation of relativistic heavy-ion A + A collisions as 
the superposition of a soft hydro-type state and a hard 
multiparton state. Both states are treated independently 
HYD JET++ [9] is a further development of its predeces- 
sors: the HYDJET [10] and FASTMC [12, 13] Monte 
Carlo generators. The soft part is based on a hydro- 
dynamical parametrization of the initial state providing 
the thermal hadronic state generated on the chemical 
(single freeze-out scenario) or thermal (thermal freeze- 



out scenario) freeze-out hypersurfaces represented by a 
parametrization of relativistic hydrodynamics with given 
freeze-out conditions [12, 13]. The mean multiplicity of 
hadron species crossing the spacelike freeze-out hypersur- 
face is calculated using effective thermal volume approx- 
imation. Note that unlike FASTMC in HYDJET++ the 
value of effective volume of the fireball V e g is generated 
for each event separately. V e e is proportional to the mean 
number of participating nucleons at the considered cen- 
trality (impact parameter b), which is calculated from the 
generalization of the Glauber multiple scattering model 
to the case of independent inelastic nucleon-nucleon col- 
lisions. In the case of the thermal freeze- out scenario the 
system expands hydrodynamically with frozen chemical 
composition, cools down and finally decays at a thermal 
freeze-out hypersurface [12]. The two- and three-body 
decays of the resonances with branching ratios are taken 
from the SHARE particle decay table [14]. 

The model for the hard multiparton part of the HY- 
D.IET++ event is the same as that for the HYDJET 
event generator. A detailed description of the physics 
framework of this model can be found in Ref. [10]. The 
approach for multiple scattering of hard partons in dense 
QCD matter, that is, quark-gluon plasma, is based on an 
accumulating energy loss, the gluon radiation, and colli- 
sional loss being associated with each parton scattering 
in the expanding quark-gluon fluid. 

The routine for generation of a single hard NN colli- 
sion, PYQUEN [10, 15], is constructed as a modification 
of the jet event obtained with the generator of hadron- 
hadron interactions, PYTHIA_6.4 [16]. The event-by- 
event simulation procedure in PYQUEN includes the 
generation of the initial parton spectra with PYTHIA 
and production vertexes at the given impact parame- 
ter, rescattering-by-rescattering simulation of the parton 
path length in a dense zone, radiative and collisional en- 
ergy loss, and final hadronization with the Lund string 
model for hard partons and in-medium emitted gluons. 
Then, the full hard part of the event includes PYQUEN 
multijets generated around its mean value according to 
the binomial distribution. The mean number of jets pro- 
duced in A + A events is a product of the number of 
binary NN subcollisions at a given impact parameter 
and the integral cross section of the hard process in JVJV 
collisions with the minimal transverse momentum trans- 
fer, p™ ln . Further details of the model can be found in 
Refs. [9, 10, 12, 13]. 



III. HOW JETS CAN DIMINISH THE ELLIPTIC 
FLOW AT LHC COMPARED WITH RHIC 

In the following we consider heavy-ion collisions at 
fixed centrality a/cr gco = [b/(2R A ) 2 ] = 42%, corre- 
sponding to impact parameter b w 1.3i?A- At this 
centrality the elliptic flow is already quite strong, and 
the total hadron multiplicity is still high enough to re- 
duce fluctuations in the high-pr region of particle spec- 
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FIG. 1: (Color online) The px dependence of total elliptic flow 
(upper row) and its hydro component (bottom row) in the 
HYDJET++ model for different hadron species. Left column: 
Au + Au collisions at yfs = 200.4 GeV. Right column: Pb + 
Pb collisions at i/s = 5.5A TeV. Centrality for both reactions 
is fixed at c = a/a gco — 42%. Data for ir° flow, shown in (a), 
are taken from Ref. [18]. 
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FIG. 2: (Color online) The px distribution for different 
hadron species in HYDJET++ calculations of Pb + Pb col- 
lisions at y/s — 5.5A TeV with centrality c = 42%. 



tra. Figure 1 shows the transverse momentum depen- 
dence of elliptic flow coefficient for the most abundant 
hadron species, namely, pions, kaons, protons, lambdas 
and sigmas, produced in 1,000,000 gold-gold collisions at 
= 200,4 GeV and in ca. 500,000 lead-lead interac- 
tions at y/s — 5. 5 A TeV. At least three features should 
be mentioned here. First, for both reactions, the ini- 
tial increase of the flow with rising px is accompanied 
by the rapid falloff at high values of the transverse mo- 
mentum. Second, the pronounced feature of the RHIC 
experimental data [17] reproduced by HYDJET++ in 
Fig. 1 is the crossing of baryon and meson branches. 
Finally, the flow at LHC is almost 50% lower than the 
flow at RHIC. We also plot recent data concerning the 
elliptic flow of neutral pions at RHIC [18] onto the model 



predictions in Fig. 1(a). The experimental data indicate 
a slight decrease of the V2 at px > 2 GeV/c and its pos- 
sible saturation at the 10% level, whereas in the HYD- 
JET++ calculations the flow in the high-pT domain is 
much weaker and does not exceed 2-3%. Therefore, the 
model predictions of the flow excitation functions V2{px) 
above pt ~ 4GeV/c should be considered as qualitative 
ones. However, what is the origin of the flow drop and 
why is the V2 at LHC significantly smaller than the RHIC 
flow? These peculiarities in the behavior of the elliptic 
flow can be explained by the interplay between the hy- 
drolike part of the spectra and the jets. 

The flow of the hydro part rises monotonically up to 
t>2 — 0.5 at pt — 6 GeV/c, while the relative contribution 
of the soft processes to the total particle multiplicity de- 
creases with rising pr, so the particles with px > 6 GeV/c 
are produced merely via hard processes, that is, jets. The 
jets themselves indicate some amount of flow due to the 
jet quenching effect. [The energy loss of the high-pr par- 
tons depends on their passing length in the anisotropic 
matter, thus giving the large yield of the high-py par- 
tons in the short in-plane (x, z) direction compared with 
that in the long out-of plane (y, z) one.] However, their 
flow is small and does not exceed 5% in HYDJET++ 
calculations even at full LHC energy. Therefore, the 
superposition of the hydro part and the jets results in 
weakening of the elliptic flow after a certain transverse 
momentum. At LHC, jets turn to dominate over the soft 
processes at much lower values of px than at RHIC, thus 
effectively decreasing the elliptic flow, t>2 HC < t^ HIC at 
Pt > 3 GeV/c. It is worth mentioning that the last result 
stems from the assumption of similarity of model param- 
eters responsible for the correct description of particle 
chemical and thermal freeze-out at both energies. (We 
have considered a naive "scaling" of the existing physi- 
cal picture of heavy-ion interactions over the two orders 
of magnitude in center-of-mass energy to the maximum 
LHC energy yfs =5.5 A TeV. The linear extrapolation 
of the model parameters in log( v / s) to the LHC can be 
found in Ref. [19].) These parameters are listed in Ta- 
ble I. Some of them will probably be modified at LHC 
more seriously compared with the estimated values. Nev- 
ertheless, even if the elliptic flow at LHC in the low-py 
region will be stronger than that at RHIC, our condition 



,LHC 



< vf mc 



will be accomplished at slightly higher 
transverse momenta, say, px _ 4GeV/c. 

Last but not least, in HYDJET++ at low px the el- 
liptic flow is strictly ordered by particle masses. Light 
particles, such as pions and kaons, have larger flow than 
heavier ones, such as protons and lambdas. On the 
other hand, the slope of the px spectra of light parti- 
cles is steeper than that of the heavy particles, as can 
be clearly seen in Fig. 2. As a result, the hydro compo- 
nent of the transverse momentum distribution of heavy 
hadrons dominates until larger values of px- For instance, 
for pions and kaons, it determines the spectrum up to 
px ~ 4GeV/c, whereas for protons the hydro part domi- 
nates until px ~ 5 GeV/c. Because of this, the fall of the 
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TABLE I: The input parameters for the HYDJET++ gener- 
ator for the two reactions in question: T c h, temperature at 
chemical freeze-out; 2th, temperature at thermal freeze-out; 
/ib, baryon chemical potential; i?™Jn, maximum transverse 
radius at thermal freeze-out; y^ng, maximum longitudinal 



flow rapidity at thermal freeze-out; j/™an, maximum trans- 
verse flow rapidity at thermal freeze-out; and ?J™ ln , minimum 
transverse momentum of parton-parton scattering. For Au + 
Au, = 200A GeV. For Pb + Pb, = 5.5A TeV. 
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TABLE II: Yields of the particles produced directly and via 
the resonance decays at the midrapidity region. Feed-down 
from weak decays of strange particles is included. For Pb + 
Pb, c = 42% and y/s = 5.5,4 TeV. For Au + Au, c = 42% 
and V« = 200A GeV. 
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FIG. 3: The pr dependence of elliptic flow of direct hadrons 
(open symbols) and of all hadrons (full symbols) produced 
in the HYDJET++ model for Pb + Pb collisions at -fH = 
5.5A TeV with centrality c = 42%: (a) protons, (b) pions, (c) 
kaons, and (d) lambdas plus sigmas. 
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elliptic flow with rising transverse momentum occurs for 
light particles at smaller px, and, at px > 4GeV/c, the 
mass ordering of the V2{pt) spectra is reversed: heaviest 
particles possess the largest flow. 

IV. THE INFLUENCE OF RESONANCE 
DECAYS 

At RHIC energies the transition to meson-rich mat- 
ter has been found. As was predicted by Hagedorn [20], 
at high energies most of the particles will be produced 
through resonance decays with shifting of the average 
mass to a heavier sector. Table II shows the contributions 
of direct and resonant production for various hadron 
species including feed-down from weak decays for the Pb 
+ Pb event sample generated with HYDJET++ at the 
top LHC energy. One can see that about 80% of pions, 
70% of protons, 60% of E and A hyperons and more than 
50% of kaons are produced from the decays of resonances. 

The effect of resonance decays should be accounted for 



FIG. 4: The same as Fig.3 but for Au + Au collisions at 
V's = 200,4 GeV. 



when one considers the development of the elliptic flow 
as well. The difference between V2 of all these hadrons 
and V2 of only directly produced hadrons is displayed in 
Figs. 3 and 4 for Pb + Pb collisions at LHC and for 
Au + Au collisions at RHIC, respectively. The degree of 
influence of resonance decays on the strength of elliptic 
flow is quite different for various hadrons. The effect is 
strong for protons and S + A hyperons, rather moderate 
for pions, and extremely small for kaons. For all particles, 
except pions, the resonances either do not alter the flow 
(kaons) or increase it for both energies in question. Note 
that in the pr region below 1 GeV/c the elliptic flow of 
direct pions is larger than the flow of all pions. This 
finding is in line with the result obtained in relativistic 
(2+1) hydrodynamics [21]. One can conclude that the 
soft pions emitted from the decays of resonances have 
lower momentum anisotropy, whereas heavier particles 
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TABLE III: Relative yields of pions and protons produced 
directly and via the resonance decays in HYDJETH — K Cal- 
culations are done for Au + Au collisions at RHIC and Pb + 
Pb collisions at LHC with centrality c=42%. For Pb + Pb, 
y/s = 5.5A TeV. For Au + Au, = 200A GeV. 
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FIG. 5: (Color online) The pr dependence of elliptic flow for 
(a) charged pions and (b) protons plus antiprotons produced 
both directly (circles) and in A decays (squares) in the HY- 
DJET++ model for Pb + Pb collisions at -*fs = 5.5,4 TeV 
with centrality c = 42%. (c) and (d) The same as (a) but 
for charged pions produced in decays of p and uj, respectively. 
The flow of resonances is shown by triangles. 



(p T aryon ) > (p™ cson ). Since the transverse momentum 
spectrum of the produced pions is softer, the pion elliptic 
flow gets an extra boost at low pt from the intermediate- 
Pt flow of heavy resonances see 

Figs. b(a)and6(a) 
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FIG. 6: (Color online) The same as Fig. 5 but for Au + Au 
collisions at ^/s — 200A GeV. 



and hard pions demonstrate larger elliptic flow of the 
resonances. 

In order to study this peculiarity in detail, let us con- 
sider the case of pion and proton flow. Figures 5 and 
6 show differences in spectra of secondary pions and 
protons coming from (anti)delta (A ++ , A + , A , A~) de- 
cay. When a heavy baryon resonance decays into a sec- 
ondary baryon plus a pion, the majority of its trans- 
verse momentum, because of the decay kinematics, is 
carried by the baryon, while the pion is produced with 
lower p T . As seen in Figs. 5(a), 5(b), 6(a), and 6(b) 
the resulting elliptic flows of baryon and pion decay frac- 
tions possess similar amplitudes (but not the same pr) 
as was carried by the resonance before the decay, that 

is, ( W ™ sonancc ) RJ (wJT ry0n ) « (i;™ son ), but (p^sonanco^ _ 



. In contrast, the secondary baryon has practically the 
same V2(pt) dependence as that of the mother particle, 
as shown in Figs. 5(b) and 6(b). 

The relative contributions of different resonance chan- 
nels to the yields of protons and pions at RHIC and LHC 
energies are presented in Table III. We see that, for ex- 
ample, lambdas and deltas alone, together with their an- 
tistates, produce about 45% of all protons and antipro- 
tons. On the other hand, the contribution of baryon res- 
onances to meson spectra is relatively weak. As follows 
from Table III, many secondary pions are produced from 
p and uj mesons. The momentum distributions of these 
pions are quite different, as can be seen in Figs. 5(c), 
5(d), 6(c), and 6(d). 

Elliptic flow of pions from the p — > tttt decay al- 
most coincides with v^, while in the three-particle decay 
uj — > 7T7r7r pions are getting an obviously softer pt distri- 
bution, and thus their elliptic flow is transferred to the 
softer pt region compared with v^^Pt)- Although the 
contributions to the elliptic flow of pions coming from 
baryon and meson resonances at pr < 1.5GeV/c effec- 
tively compensate each other, the resulting elliptic flow 
appears to be a bit lower than that of directly produced 
pions. For higher transverse momenta the contribution 
from the decays of K mesons and heavy resonances de- 
termines the observed excess of the resulting pion flow 
over the direct pion flow. 

In general, the resonance contributions sometimes in- 
crease and sometimes decrease the initial elliptic flow 
assigned to directly produced particles. Some hadrons, 
such as the <f> meson, do not get a feed-down from the res- 
onances, and their flow profiles remain unchanged. Mod- 
ification of the elliptic flow of pions and light baryons, 
especially pronounced at LHC energy, can lead to viola- 
tion of the hydro- induced mass hierarchy in the V2(pt) 
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FIG. 7: (Color online) (a) The pr dependence of elliptic flow 
in the HYDJET++ model for different hadron species (lines) 
and comparison with RHIC data (symbols), (b) The same as 
(a) but for scaling variables, Vijnq vs KEt- 



sector. 



V. NUMBER-OF-CONSTITUENT-QUARK 
SCALING 

One of the most interesting features observed in the de- 
velopment of hadron elliptic flow at RHIC is the so-called 
number-of-constituent-quark (NCQ) scaling [22, 23]. It 
appears that elliptic flow of any hadron species depends 
on the transverse kinetic energy KEt = tut — mo in a 
similar manner provided both V2 and KEt are divided by 
the number of constituent quarks, n q , in a given hadron; 
thas is, n q — 2 for a meson and n q — 3 for a baryon. 
The observance of the NCQ scaling in a broad kinematic 
range implies the formation of elliptic flow already on 
a partonic level. Recent experimental studies based on 
higher statistics indicate that the scaling holds up only 
until KE T /n q fa 1 GeV [17]. 

The experimental data on particle elliptic flow and 
NCQ scaling in 1,000,000 minimum bias Au + Au col- 
lisions at ^/s = 200^4 GeV are shown in Fig. 7 together 
with the results of the HYDJET++ simulations. Since 
no information about quark content is used in the hydro 
part of the HYDJET++ model, the scaling might not 
necessarily have been observed. For instance, a <p meson 
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FIG. 8: (Color online) Upper row: The KEr/riq dependence 
of elliptic flow for (a) direct hadrons and (b) all hadrons pro- 
duced in the HYDJET++ model for Au + Au collisions at 
y/s = 200 A GeV with centrality c = 42%. Bottom row: The 
KEr/riq dependence of the ratios (v2/n q ) j (v% /3) for (c) di- 
rect hadrons and (d) all hadrons. 



is heavier than a proton and, according to the mass or- 
dering, its normalized flow v^pt)/^ is larger than the 
proton flow v^ipr)/^ already at pt > 0.8 GeV/c. On the 
other hand, the elliptic flow in the model is fit to describe 
experimental data in the low-p^ range, as displayed in 
Fig. 7(a). Therefore, the dependence vijn q (KET /n q ) 
may also exhibit the scaling trend. The model calcula- 
tions presented in Fig. 7(b) seem to obey the approximate 
NCQ scaling up to KET/n q as 0.7 GeV or even to higher 
values if we exclude heavy hyperons, such as 2 and S, 
from our consideration. (There are several experimental 
and theoretical lines of evidence that multistrange hy- 
perons are frozen earlier; see Ref. [24].) Is this just a 
coincidence, and what is the role of resonance decays? 

To answer these questions we plot HYDJET++ exci- 
tation functions V2/n q (KET /n q ) for direct hadrons only 
and for all hadrons produced in gold-gold collisions at 
top RHIC energy in Figs. 8(a) and 8(b), respectively. 
In order to see possible deviations from the scaling be- 
havior more distinctly, the ratios v% /n q : w^/3, i = 
TT ± ,K ± ,p/p,A/A,(j) are presented in Figs. 8(c) and 8(d). 
One may note that, for direct hadrons, the NCQ scaling 
is fulfilled within the 20% accuracy limit in the interval 
0.2 < KEt < 0.8 GeV. At higher and lower transverse 
energies the curves diverge, as shown in Figs. 8(a) and 
8(c). The situation is markedly improved after the res- 
onance decays are taken into account, as demonstrated 
in Figs. 8(b) and 8(d). All ratios, except for kaons at 
KEt < 0.2 GeV, are much closer to unity, and the scal- 
ing appears to hold at 10% up to KE T > 1.0 GeV. The 
spectrum and, therefore, the elliptic flow of <f> mesons 
remain unchanged because of the absence of resonance 
feed-down. In contrast, protons and lambdas are signifi- 
cantly boosted by heavy resonances, and their flows are 
enhanced at KE T > 0.5 GeV. 



7 



"~>J 0.06 



(a) direct hadrons 


(b) all hadrons 


- 




— i — i — i — i — i — i — i — i — i — i — i — f t 

L(c) / / 


- 1 1 1 1 i 1 1 1 1 i 1 1 1 1 f 

\ (d) 'K ; - J 











% 0.5 1 1.5 0.5 1 1.5 

KEp/riq [GeV] 



FIG. 9: (Color online) The same as Fig.8 but for Pb + Pb 
collisions at y/s = 5.5^4 TeV. 



The same mechanism also works at LHC energies. 
However, because of the influence of jets at relatively 
low transverse momenta, there is no resemblance of the 
NCQ scaling for directly produced hadrons, as shown in 
Figs. 9(a) and 9(c). Here, all curves seem to diverge in the 
entire transverse energy range. Again, the decays of reso- 
nances essentially raise up the elliptic flow of protons and 
lambdas displayed in Figs. 9(b) and 9(d). However, the 
effect of jets at LHC is too strong and, although all ratios 
of particle flows presented in Fig. 9(d) are almost paral- 
lel within the energy range 0.1 < KEt < 1 GeV, the re- 
alization of theapproximate number-of-constituent-quark 
scaling becomes worse compared with the RHIC case. 

A few remarks on results obtained in other calculations 
are in order here. The effect of resonance decays on the 
elliptic flow of stable hadrons in Au + Au collisions at 
RHIC was studied in Ref. [25] in the framework of the co- 
alescence model and in Ref. [26] within the parametrized 
thermal model. In contrast to our study, the authors 
of both works postulate that the elliptic flows of di- 
rectly produced stable particles and resonances obey the 
number-of-constituent-quark scaling. After that, the res- 
onances account for a substantial part of deviation of pion 
elliptic flow from the NCQ scaling. We do not assume 
the automatic performance of the NCQ scaling, but in 
our model calculations the resonance decays also modify 
the pion elliptic flow by enhancing it in the intermediate- 
Pt region. In accordance with Ref. [25] kaon elliptic flow 
in our calculations is not affected by the decays of res- 
onances, but we see the strong contribution from heav- 
ier resonances to the momentum anisotropies of protons 
and lambdas. This outlines the importance of a com- 
plete table of resonances in any microscopic or macro- 
scopic model designed for the description of heavy-ion 
collisions. For instance, the present version of the HYD- 
JET++ model employs ca. 350 baryon and meson states 
as well as their antistates. 

It was suggested recently [27] that NCQ scaling could 
be used as a unique probe of strongly interacting partonic 



matter. The scaling holds if the hadrons are produced 
predominantly from the coalescence of quarks, whereas it 
is broken if the hadronization proceeds via the string frag- 
mentation. We see, however, that the jet effects become 
increasingly important at ultrarelativistic energies. Par- 
ticularly, the interplay between jets and hydrolike flow 
can cause breaking of the NCQ scaling for the particle 
elliptic flow at LHC. 



VI. CONCLUSIONS 

The elliptic flow pattern in Pb + Pb collisions at y/s = 
5.5A TeV and in Au + Au collisions at y/s = 200A GeV 
is analyzed for different hadron species in the framework 
of the HYDJET++ Monte Carlo model. The model con- 
tains both parametrized hydrodynamics and jets. This 
allows one to study the interplay between hard and soft 
processes and to reveal their role in the formation and 
evolution of the elliptic flow. It is worth mentioning that 
the model parameters for the soft hydro component are 
fixed at RHIC energy to simultaneously describe hadron 
yields, energy spectra, anisotropic flow, and femtoscopic 
momentum correlations. These parameters are not ex- 
pected to change significantly at the LHC energy. 

Three general tendencies have been observed. First, 
at a certain transverse momentum, particles produced 
via jets start to dominate over the hydrolike particles. 
Since the flow of jet particles determined merely by the 
jet quenching is weak, the resulting flow V2(pt) expe- 
riences falloff in the intermediate pt range. Abundant 
production of jets at LHC will effectively decrease the 
hadron elliptic flow at smaller values of px compared 
with RHIC. Therefore, we expect that the elliptic flow 
of stable hadrons at LHC will be smaller than its RHIC 
counterpart, v^ HC (p T ) < vf HIC {p T ), at p T > 3GeV/c. 

Second, jets also account for changing of the mass or- 
dering of the hadron elliptic flow at intermediate and high 
transverse momenta. After the falloff the mass ordering 
of the V2{pt) spectra becomes opposite to the initial one; 
that is, the heavier particle has the larger flow. 

Third, resonances are shown to significantly modify the 
elliptic flows of directly produced hadrons. The ^(fr) of 
kaons is not affected by the decays, while the flows of pi- 
ons and especially protons and lambdas are enhanced. It 
is the effect of resonance decays that pushes the particle 
spectra toward the fulfillment of the constituent quark 
number scaling. At the LHC energy the jets should alter 
the hadronic elliptic flow already at intermediate trans- 
verse momenta, thus completely removing resemblance 
of the NCQ scaling for directly produced hadrons. After 
the resonance decays the elliptic flows of different hadrons 
move closer to each other; however, only the approximate 
KEr/nq scaling holds. 
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